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INTRODUCTION 


This paper examines the relationships between several semiempirical 
fracture analyses and the R-curve concept of fracture mechanics. These 
relationships may explain why a semiempirical fracture analysis will yield 
good results with one set of data and poor results with another. They may 
also indicate which analyses deserve further consideration. 

Over the past decade a number of semiempirical fracture analyses have 
been presented (refs. 1-5). These analyses all attempt to correlate failure 
stresses for precracked tension specimens with initial crack length over a 
range of crack lengths. The correlations involve the determination of one 
(refs. 1-3) or two (refs. li,5) empirical parameters from test data. The 
parameters are treated as material properties which are independent of the 
specimen and crack configuration but which are functions of specimen thick- 
ness and such variables as heat treatment and test temoerature. The analy- 
ses do not always provide good correlations using data sets other than those 
chosen by the original authors. To date these analyses have only been 
formulated for and applied to test specimen configurations. Thus their 
applicability to the design of complex structural configurations is uncertain 

The progressive development of the R-curve concept has been reviewed 
in ref. 6. The concept postulates that, for a given material and thick- 
ness, there is a unique relationship between the amount of stable crack 
growth under rising load and the crack-tip stress intensity factor. This 
relationship is called the crack-extension resistance curve, or R-curve, and 
represents the response of the material in the vicinity of the crack tip to 
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externally imposed loading. If the R-curve is known, both failure load and 
critical crack length can be predicted (as functions of initial crack length) 
for any specimen or structural configuration for which an appropriate stress 
intensity analysis is available. Thus the R-curve concept appears to be a 
more useful method than any of the semiempirical analyses. 

If the R-curve for a given material and thickness is available, one can 
calculate fracture stress as a function of original crack length for any 
test specimen configuration. The converse should also be true. That is, if 
a relationship between fracture stress and original crack length is avail- 
able, one should be able to calculate the corresponding R-curve. This 
observation was the impetus for the present study, which was undertaken to 
test the following hypotheses: 

1) For each semiempirical fracture analysis (SEFA) there is an equiva- 
lent R-curve (ERC) whose magnitude and shape are determined by the SEFA 
formulation and its empirical parameters. The ERC is equivalent in that it 
predicts exactly the same relationship between fracture stress and original 
crack length as the SEFA. 

2) A SEFA will correlate residual strength data (fracture stress against 
original crack length) closely if its ERC closely matches the actual R-curve 
of the material in question, and will correlate poorly if the match is poor. 

This paper first reviews some characteristics of the R-curve concept 
when applied to finite-width specimens. Next the conditions for equivalence 
between a semiempirical analysis and an R-curve are derived. A hypothetical 
material is emoloyed to study the relationship between R-curves and semiem- 
pirical analyses. Finally, equivalent R-curves are developed for real 
materials using data from the literature. 
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SYMBOLS 

a Length of single-tip crack or half-length of double-tip crack, 

equals a 0 + A 

n 

E' Effective modulus, ec;uals E for plane stress or E/Cl-^' 1 } for 

plane strain, where E is Young's modules and k is Poisson's ratio 

C-^ Strain energy release rate 

G- Crack extension resistance 

G c Strain energy release rate or crack extension resistance at instati lit 

Kj Opening-mode stress intensity factor 

r. lumber of crack tips (one or two) 

'■J Specimen vddth 

Y Stress intensity calibration factor, Kj/CT/a" , a dimensionless 

function of ,A 

(X Sensitivity factor, eq. (1) 

A Effective crack extension (sun of ohysi.cal crack extension plus 

a plastic zone correction) 

A Relative crack length, n&A r 

CT Stress normal to crack 

6/^ Ultimate tensile strength 

>7y S Yield strength 

Subscripts: 

c at critical or instability condition 

o original value prior to loading 

R -CHIVE CONCEPT 

The R-curve concept and its historical development are reviewed in ref. 

6. A general representation of the concept is shown in figure la. The 
strain energy release rate is given b"- the expression 
G a = Y 2 (T 2 a/E' 

and represents the driving force (per unit thickness) tending to cause crack 
propagation. The material's resistance to crack propagation, G, , is a 
function of crack extension, A . At the critical stress 0~ c the driving 
force curve and the R-curve are tangent. Beyond the point of tangency the 
driving force increases faster with crack length than does the material's 
resistance. This instability condition represents the failure of the body. 
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The point of tangency defines the fracture toughness, G c , and the critical 
crack length, 2a c . For an infinite body, Y is a constant and the driving 
force curve is a straight line Thus both the fracture toughness, G c , and 
the amount of crack extension at instability, a c > increase vrith increasing 
original crack length. 


For cracks in simple finite bodies and test specimens, the trends are 
somewhat different. If we define a dimensionless sensitivity factor as 


o( * 


A dY 
Y d\ 


( 1 ) 


then the crack driving force curve and its slope (for constant stress) are 
E'G a * Y 2 <T 2 a 

E'“ = Y 2 </ 2 (1 + 2c< ) 

For convenience, the crack extension resistance curve and its slope are 
written here as 


g(A) B E'Gp_ 
g ' ( A ) 5 E' dG p /dA 


At the instability point, G A =G^ and dG A /da=dC-jj/dA (see figure la). I i 
g( A ) and g 1 ( A ) are mathematically describable, the instability point 
is determined by the simultaneous solution of two equations, 


F'G, 


= v 2 
"C 


°~C. ( a o + A c ) = g( A c ) 


,11 

da 


E c(l + 20 <c) - S’(A C ) 


( 2 ) 

(3) 


The coefficients Y and 0( are usually expressed as trigonometric or 
polynomial functions of the relative crack length A . As a result, a 
closed-form simultaneous solution is seldom possible, and numerical 
methods must be used to solve for A c . Then G c = g( A c )/E' and the 
fracture stress <T C is determined from eq. (2). 
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As the initial crack length is increased from zero, both G c and A c 
increase. However , due to the fact that dl/dA continually increases with 
>. , both G c and A c reach maximum values which depend on the specimen 
width, ™ , and the forms of both the driving force curve and the P.-curve, 
As a Q is increased still further, both G c and A c begin to decrease. 
This behavior is shown schematically in figure lb where instability curves 
are shown for a wide range of initial crack lengths. The locus of all 
instability points is shown by the dotted line. From figure lb it is also 
apparent that there are pairs of original crack lengths, say (a 0 )i « 0.2 W 
and (a o )j«0.? W , which will have the same critical crack extension, 
(A c )l,2 > and fracture toughness, (° c )i^2 * From oration (2), the 

fracture stresses for these original crack lengths are related by 

(Qg (Ych / ( a o) 1 + ^c)l,2 

(^c)l " ( y c) 2 » ( a o) 2 + ( A c)l,2 

Thus there is a relationship between fracture stresses for short cracks and 
long cracks which is implicit in the R-curve concept, and this relationship 
is a function of the specimen type and the shape of the R-curve. 

It should be noted that, in this paper, A is the effective crack 
extension. It is the sum of the physical crack extension plus an adjustment 
to account for the effect of crack- tip plasticity. 

EQUIVALENCY ANALYSIS 

In the preceding section it was shown that, if a mathematical formu- 
lation of the R-curve is available, fracture stress can be determined as a 
function of original crack length. In this section it will be shown that, 






if an equation for fracture stress as a function of original crack length 
is available, the equivalent 7 .~curve can be determined. In this paper, 
only the case where specimen width V is constant is considered. Parallel 
derivations for the cases where a 0 is constant and where & 0 A'l is 
constant are given ir. ref. 7 3 which contains a more complete treatment of 
this same subject. 


Differentiating both sides of eq, (2) with respect to A c , 


* i 0 2 <r c 2 -2« 0 [^*i] * Y c s (vA:>a§- 0 ( <r c > 2 + 
tl «>M^ +l ] ' g,(A c) 


and substituting eq. (3), results in 

0 . <r 0 2 [(X+2fi< 0 ) ♦ (v/g g§- c «T C 2 ) 

Assume that there is a function f such that we can define 
f(a 0 ) e ^T c 2 


f'(a 0 ) ■ d( (f c 2 )/da Q 


for W * constant 


Then eq. (I4) becomes 


0 « 


da 0 

dA. 


[(1 + 2C t c ) f(a 0 ) + (a 0 + A c ) f'(a 0 )] 


W) 


and since da 0 /dA c / 0 we have 


(1 + 2tf c ) f(a 0 ) + (a 0 + A 0 ) f'(a 0 ) 


0 » 


( 5 ) 
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For cracks in infinite bodies, 0( •= 0 


A c 


f(a 0 ) 
-f 1 (a 0 ) 


and eq. (?) 


V A }}! 

bJL?°° R 


(6) 


which, after substituting the infinite-body formulations for f(a 0 ) and 
f (a 0 ) from the APPENDIX, gives A c for any a Q in terms of the empirical 
parameters. Then terms can be rearranged to give a 0 as a function of A c , 


sqy 

a Q - F( A c ) 

Substituting this function into eq. (2) yields 


E'G C = Y c 2 [f( A c ) + A c ] »f[F(A c )] (7) 

Since eq. (6) gives A c for any value of a 0 , eq. (7) must give E ! G C for 
any and all values of A c , which is a definition of the R-curve. Thus, 
after writing the function F in terms of the empirical parameters, it is 
aporopriate to write eq. (7) in the general terms of E’Gr and A , rather 
than E'G C and A c . The end result is an explicit ERC formulation in 
terms of the empirical parameters. 

To determine the ERC for cracks in finite bodies, an indirect method is 
required. First, the finite-body formulations for f(a 0 ) and f'(a 0 ) from 
the APPENDIX are substituted into eq. (?), Because of the more complicated 
nature of the finite-body formulations, it is unlikely that an explicit 
function F(A 0 ) will be obtainable. But for any given value of A c , a 0 
is a root of eq. (?) which may be found by standard numerical methods and 
which represents a single value of f( A c ) . Substituting this value into 
eq. (7) yields a discrete point on the ERC. By incrementing A c and repeat- 
ing the calculation, the ERC can be determined point by point. 
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Dimensionless Equivalent R-Curves 


It is helpful and more efficient to first compare semiempirical fracture 
analyses (ST T A) and equivalent R-curves (ERG) on an analytical basis. This 
is most easily done using the oroblem of a crack in an infinite olate as a 
baseline . 

The infinite-body ERC for Kuhn's analysis (ref. 1) is obtained using the 
method described in the paragraph containing eqns. (6) and (7), Substituting 
eons. (Al) and (A2) from the A^PSmx into eqns. (6) and (7) with Y c ^= TT 
yields 



(3a) 


E'G R 


rrcr u 2 a 

1 * C m 2 A 


(9a) 


where C m is an empirical parameter having units (L - ^). Equation (9a) is 
plotted in dimensionless form in figure 2a. This curve obviously resembles 
an R-curve and might be expected to closely match some (but not all) 
experimental R-curves . 

The infinite-body ERC for Orange's analysis (ref. 2) is obtained in the 
same manner. Using eqns. (A3) and (Ah) results in 

A c » {Y^/(T u ) 2 /7r (8b) 

E'G h = K u 2 (9b) 

where Ky is an empirical parameter having units (FL“^^ 2 ). These equations 
define a single point. In order to relate this single point to the R-curve 
concept, the point may be thought of as the corner of a step-function, and 
that step-function might in turn be considered as a very simple approximation 


of an actual R-curve. 
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The infinite-body ERC for Feddersen’ s analysis (ref. 3) is obtained by 
substituting eqns. (a£) and (A6) into eqns, (6) and (7) as before, yielding 



2 

for a n ! ( 9/U?r)( K /(T. b ) where K„ is an empirical parameter having units 
^ c ys c 

(FIT 3/2 ). Equation (8c) requires that critical crack extension decrease as 
original crack length increases from zero. This is in direct opposition to 
the R-curve concept and is not supported by any data known to this author. 
Equation (9c), which is plotted in figure 2b, does not look at all like an 
R-curve but does satisfy the requirements of coincidence and tangency. For 
a Q =0 , the point of tangency is the right-hand terminus of the curve. As 
a 0 increases, the point of tangency moves downward and leftward along the 
curve. Finally, at a 0 =( 9/U-R')(K c /<^r S ) , the point of tangency is the 
left-hand terminus. 


The infinite-body ERG for Newman's analysis (ref. lj) is obtained by 
substituting eqns. (All) and (A12) into eqns. (6) and (7) as before, yielding 



n K f , 

1 a 0 


^c 

= tf-u V 

IT 

(8d) 

E'Gr 


CT u 2 A 

(9d) 

m^/f 7" 

+ «r u Ar) 2 A 


where is an empirical parameter having units (FL“3/2) and m is a 

dimensionless empirical coefficient which is not greater than unity. Equa- 
tion (9d) is plotted in dimensionless form in figure 2c. This ERG is 

p 

asymptotic to E'Gpj^K^ , and the coefficient m determines the rapidity of the 
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approach. As m decreases from unity to near-zero, the ERC develops a 
progressively sharper knee. The flexibility of this two-parameter KRC 
should allow it to match R-curves for a wide range of real materials. 


The infinite-body ERG for Bockrath's analysis (ref. £) is obtained by 
substituting eons. (A13) and (Allj) into eqns. (6) and (7) as before, 
resulting in 

A c »■ (oj/2)a 0 (8e) 

m 

E'Gr s ^(2 + ( — A (9°) 

where to is a dimensionless empirical coefficient and Kj c is an empirical 
parameter having irrational units (FI.**). Equation (9c) is plotted in figure 
2d. This SRC has no asymptote, and its slope is infinite at A>=0 . Except 
for notation, it is identical to the R-curve model proposed by Broek as 
equation (10) of ref. (8). Broek' s model was derived using R-curve concepts 
and the experimental observation that, for small cracks in wide soecimens, the 
critical crack length is often proportional to the initial crack length. 


At this point ve can state the following. For each SEFA, in its 
infinite-body form at lerst, there is indeed an ERC. For four of the 
analyses considered, the ERC resembles or approximates an actual R-curve. 
The ERC for Feddersen's analysis does not resemble an R-curve and will not 
be considered further. 


Comparisons Using Synthetic Data 


Hypothesis (2) of the INTRODUCTION postulates that a SEFA will 
correlate residual strength data closely if its ERC closely matches the 
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actual R-curve and will correlate poorly if the match is poor. To test this 
hypothesis we would need, as a minimum, both residual strength data (for 
several specimen sizes having a wide range of crack lengths) and R-curve 
data for two materials having significantly different R-curve shapes. Since 
no such body of data is known to this author, it was necessary to synthesize 
one. 'This was done by formulating two R-curve equations using the following 
guidelines. lirst, to avoid exact fits, neither equation should be mathema- 
tically equivalent to one of the SRC formulations previous^ derived. Second, 
one R-curve should have a definite, knee, the other should be gently curving. 
Using those equations, synthetic test data can be generated by instability 
analysis for any specimen size and typo. An advantage of this approach is 
the total absence of data scatter. 


Unobtainium is assumed to be a heat-treatable material. Since the 
material is imaginary, the units will be left to the reader's imagination. 
In the annealed condition, its ultimate tensile strength is 1 £0 and its 
R-curve is given by 

E'G r = 8,000 \flOA - A 2 (10a) 


In the aged condition, its ultimate tensile strength is 200 and its R-curve 
is given by 


£ 0,000 , . . 

E'Gft - — p — arctan (10A) (10b) 
These are shown in figure 3- The coefficients in eqns, (10) were selected 
so that the significant features of both curves would lie within the ranges 
0< 2£,000 and 0£A*1. 


The pseudotest data points are calculated using conventional instability 

analysis as follows. Dividing eq. (2) by eq. (3) and rearranging terms gives 

s(A c ) a o + ^c 
s r (A c ) “ 1 '+ 2C<, 


0 


( 11 ) 
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The functions g(A c ) and g'(A c ) are given by one of eqns. (10) and its 
derivative. The factor o( c is determined using eo_. (l) and the secant 
stress intensity calibration factor for ur.iformly-loaded center-crack speci- 
mens (ref. ?). Then, for prescribed values of a 0 and W , A c is the least 
positive root of eq. (11). This root can be found by any of several numeri- 
cal methods. Next, A c is substituted back into eq. (10) to calculate E'0 C 
Finally, the fracture stress <T is obtained from eq. (2). 

The "specimens 11 that are studied here were sized as follows. For the 
infinite-width pseudotests, the initial crack half-lengths a 0 were chosen 
(by trial and error) to give A c values well distributed over the entire 
ll-curve. For the finite-width pseudotests, the specimen widths were fixed 
at 6 times (first series) and U times (second series) the largest initial 
crack half-length in the infinite -width series. The calculated values of 
stress and crack extension at instability are given in Table I. The values 
of stress and initial crack length were then used as inputs to the various 
semi.empirical analyses. 

The empirical parameters were determined as follows. Kuhn's parameter 
C m was calculated for each specimen in the infinite-width series using eq. 
(Al). The sinole average of seven values, C m , is given in Table I. The 
bar is used here to denote the average value for one data set. Orange's 
parameter, K u , was also calculated for each specimen in the infinite-width 
series using eq. (A3). The average value given in Table I is a weighted 

average determined in the same manner as eq. (6) of ref. 2. Newman's para- 
meters and m were determined using the least-squares procedure given 

in Appendix C of ref. ii. Bockrath's parameters Krp c and to were deter- 
mined by a least-squares fit of eq. (Al3). Since Bockrath's method is 
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restricted to cases where the crack area is less than 10 percent of the 
gross area, specimens having a o c>W/20 were excluded from the least- 
squares fit.. 

The equivalent R-curves were calculated as follows, lor the infinite- 
width series, the empirical parameters from Table I were simply substituted 
into the aopropriate one of eqns. (9). For the finite-width series, the 
indirect method (described in the paragraph containing eq, (7)) was used. 
Specifically, eqns. (A 9 ) and (AlO) or eqns. (A13) and (Alii) were used along 
with eqns. ($) and (7) and the appropriate empirical parameters. 

The residual strength of the infinite-width Unobtainium is shown in 
figure h. For the annealed condition, Bockrath's semiempirical fracture 
analysis (SEFA) provides a nearly perfect fit to the pseudodata. When 
ranked according to the sum of the squares of deviations, Newman's SEFA, 
Kuhn’s, and Orange's follow in that order. For the aged condition, the 
ranking is quite different. Fere Newman's SEFA provides a nearly perfect 
fit, with Orange's, Kuhn's, and Bockrath's following in that order. The 
equivalent R-curves (FRO) are shown in figure ?. For the annealed condition. 
The Bockrath SRC matches the actual R-curve almost perfectly. When ranked, 
according to the integral of the sauare of the deviation, the Newman ERC, 
the Kuhn ERC, and the Orange ERC follow in that order. For the aged condi- 
tion, the Newman ERC is the best match to the actual R-curve. The Bockrath 
SRC, the Kuhn ERC, and the Orange ERC follow in that order. The Orange ERC, 
although crude, is a better approximation for the aged material than for 
the annealed. 

Since the Orange ERC is rather crude and since Kuhn's SEFA is equivalent 
to a special case of Newman's (see APPENDIX, following eq. (A12)), these two 
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were not considered further. The Newman and Bockrath)EN©OCffii : 'Qf0AlSI?5j'* : ' e “ 
width center-crack series are shown in figures 6 and Hera the same 
trends are seen as in the infinite-width series. The Bockrath ERC is 
the better match for the annealed condition, while the Newman ERC is the 
better match for the aged condition. Note in Table I that the empirical 
parameters K Tc , <Z> , K f , and m all vary slightly with specimen width. 
The EHCs shown in figures £-7 are also distinctly different for different 
specimen widths, but the differences are slight. 

The results of this exercise using synthetic data can be summarized as 
follows. Hypothesis (2) of the INTRODUCTION postulates that an SEPA will 
correlate residual strength data closely if its ERC closely matches the 
actual R- curve and will correlate poorly if the match is poor. Strictly 
speaking, this hypothesis cannot be proven, since the ERC magnitude and 
shape depend on empirical parameters which must be obtained from residual 
strength data. However, the converse appears to be true. That is, if a 
SEFA correlates residual strength data closely, its ERC will closely match 
the actual R-curve. Furthermore it is apparent that if, for a given 
material and thickness, the R-curve is unique the various empirical para- 
meters are not, and vice versa. 

COMPARISONS USING ACTUAL TEST DATA 


As mentioned earlier, experimental studies containing both residual 
strength and actual R-curvo data are relatively few in number. Nevertheless, 
enough were found in the literature to allow some comparisons to be made 
using actual data obtained from real materials. 


NASA Data for 2011 i-t 6 Aluminum Alloy 


In reference 10 this author presented tost data for 2011 i-T 6 aluminum 
alloy specimens 1,? mm (0.06 in) thick, tested at 77 K (-320°F). 'figure llj 
of that reference presented typical curves of crack growth against applied 
stress for notches of six initial lengths in 30 -cm (12-in) wide specimens. 
Those curves were developed by slotting individual crack growth data points 
for replicate specimens, then drawing a smooth curve to give a good visual 

average. For the present report, those data were re-analyzed. The crack 
extension resistance and effective crack length were computed for each data 
point as 

E'G r = <T 2 7Ta e ££ secant(57a ej [-£/W) 
a e ff = a + (B'Gr/2 1t<r£) 

respectively. Since these equations are transcendental, an iterative solu- 
tion was required. A total of 176 data points were obtained from 17 speci- 
mens with initial crack lengths (2a 0 ) ranging from 3 mm to 100 mm (1/8 to 
li in) . The empirical parameters for the Newman and Bockrath SEFAs were 
determined in the manner described earlier. As before, only specimens with 
a 0 < V//20 were included in the Bockrath analysis. The fitted empirical 
parameters are listed in Table II, 

Residual strength is shown in figure 8a. Newman’s SEFA gives a good 
fit over the entire range, Bockrath 1 s SEFA fits the short-crack, data fairly 
well, but the fit would be poor if extrapolated to longer cracks. The 
R-curve data points and the equivalent R-curves are shown in figure Ob. 

Both ERCs fit the data rather well, with Bockrath 1 s somewhat better at small 
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crack extensions and Newman’s somewhat better for larger extensions. These 
results suggest that the ERG concept applies to real data as well as to 
synthesized data. 


PAGE ^ 
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Boeing Data for 2219-T8? Aluminum Alloy 

Earlier it was shown using synthetic data that, if the actual R-curve 
is unique, one obtains slightly different values of the empirical parameters 
from data sets for specimens having different widths. Data from the Boeing 
Co, for 2219-T87 aluminum alloy specimens have the same characteristics. 
These data originally appeared in ar, internal report (Eichenberger, T. W. : 
Eracturc- Resistance Data Summary, Report D2-209i;7* Boeing Airplane Go., 

June 1962), but are also tabulated in reference 1. Center-crack specimens 
2.9 mm (0.10 in) thick were tested. The data for specimens 60 and 120 cm 
(?h and hf in) wide are used here because they cover a ”ide range of 
initial crack lengths. 

Bockrath's analysis was not applied since only one of the wider speci- 
mens had a. o <W20 . Newman's parameters were determined separately for 
each specimen width, and somewhat different values were obtained as can be 
seen in Table IT. The residual strength curves fit the data quite well, as 
can be 3i»cn in figure 9a, with the average error being less than percent. 
Using a method that is outside the scope of this paper, it was found that 
the. actual R-curve for this material could be estimated by 

E'G r = a.OTxlO 1 ^ A 0 *^ 

where E'G R is in N^/ui^ and A is in cm , or by 
E'G r = 11.2xl0 9 A 0 *^ 
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where E'O^ is in lbvirr and A is in inches. Residual strengths calculated 
from this equation using conventional instability analysis are also shown in 
figure ?a. The agreement is very slightly better than for Newman's SEFA, the 
average error being less than 3 percent. The estimated R-curve and the 
Newman ERCs are shown in figure 9b, and as expected the differences are small. 
Although these data tend to support the concept of a unique R-curve, the 
difforences are so small as to be within the bounds of probable data scatter. 

CONCLUDING REMARKS 

The results of this study lead to the following conclusions: 

1. - For each semiempirical fracture analysis (SEFA) there is an equi- 
valent R-curve (ERC) whose magnitude and shape are determined by the SEFA 
formulation and its empirical parameters. The ERC is equivalent in that it 
nredicts exactly the same relationship between fracture stress and initial 
crack length (residual strength) as the SEFA. 

2. - If, for a given set of data, a SEFA correlates residual strength 
closely, its ERC will closely approximate the effective R-curve of the 
test material. 

3. - Of the five SEFAs examined, Newman's (ref. li) appears to be the 
most generally useful. Bockratb's SEFA (ref. $) , which is only formulated 
for quasi-infinite bodies, is too restrictive for widespread use. Three 
(refs. 1-3) do not appear to warrant further consideration. 

Li. - If the effective R-curve is indeed unique, then the various 
empirical parameters cannot be constant, and vice versa. 
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The analytical comparisons made herein indicate that the variations in 
Newman's parameters are small enough that the differences may well be within 
the range of normal data scatter for real materials. Thus a very carefully 
planned and conducted experiment would be required to determine which concept 
(R-curvo or SEFA) is more appropriate. 
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APPENDIX 

SEMIEMPJRICAL FRACTURE ANALYSES 


p . Kuhn (1968). - Equations (3) and (h) of ref. 1 give the fracture 
stress for a finite-width center-crack plate. For an infinite plate, those 
may be rewritten and differentiated as 

f ( a o ) c ^u 2 [l + Cm (Al) 

I’ , ( a o) = -i‘( a o) + C m /a Q ) ^ * Cj//^ (A2) 

where C m is an empirical parameter having units ( L~^) , 


T. Orange (1969). - Equation (0) of ref, 2 gives the fracture stress 
for a finite-width center-crack oLato. For an infinite plate, this reduces to 

* < a 0 ) B K u 2 jVa 0 + (K U /CT U ) 2 ] - 1 ( A3 ) 

f'(a 0 ) - -f(a 0 ) [a 0 + K u 2 /^^u 2 ] -1 (Ah) 

where K u is an empirical fracture toughness parameter having units (FL - ^ 2 ). 


C. Feddersen (1970). - For an infinite plate, cons. (6) and (10) of 
ref. 3 reduce to 



(A5) 

(A6) 


for a 0 <(9/M7)( K c /<T y g) 2 and «q. (7) to 

f(a 0 ) = K c 2 /iTa 0 
f( a o) = -f(a o )/a 0 


(A7) 

(AS) 
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Kuhn (19<»8). - Equations (3) and (L) of ref. 1 give the fracture 
s*r»sr for n flnl te-vidth center-crock plate. "or an infinite olate, those 
-.ay be rewritten and differentiated as 

f( f 1 * c n «^]‘ 9 (Al) 

f(a 0 ) - -f(a c ) * (1 ♦ C n {T 0 )' 1 ‘ C (A2) 

Ju 

vrhere C n is an empirical oarameter having units (L ) . 

T. Qrar.c» (196°). - Equation (S) of ref. 2 gives the fracture stress 
for a finite-width center-crack olnte. ? or an infinite olate, this reduces to 

f' »o> • K u 2 [*»c ♦ (K U /<T U )- ] “ X (A3) 

f(a 0 ) - -f(«o)[« 0 ♦ KbV^u 2 ]’ 1 (Ali) 

where K u is an emoirical fracture toughness parameter having units (FL“^~). 


C. Feddersen 

(1970). 

- For 

an infinite 

olnte, ecns. 

(6) and (1C) of 

3 reduce to 







f (»„ ) 

• V 

1 - 

UTT / 

27 \ K c / 

•o 

2 

(A$) 

f(a 0 ) 

■ -f(*o 

) f”— 

£17 

.(JStf. 

'Ws/ 

ao 

2 _ 

-1 

(A6) 

a 0 <(9 A*0( * 


and 

eq. (7) 

to 



f(a 0 ) 

- K c ? /*Ta 0 




(A7 ) 

f(»o> 

■ -f(«o 

)/a 0 




(A3) 
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for » 0 >(9/u^r)' K c /<Ty a ' , where <Ty a is the material's yield s*rength 
and K, is an emoirical fracture toughness parameter . 


J. Vevnan (1972). - Equation (1?) of ref. L for a finite-width center- 
crack nlate can be rewritten and di fferent iated as 


f( a 0 ) • K. 


_ V f 

*0 ■*«(*'«0 / V ) ♦ r - -*r 

1- A 0 W U 


-2 


(A9) 


f(a 0 ) - - 



which, for nn infinite plate, reduce to 

f( a 0 ) * K f 2 [/^ ♦ mK f /<7 u ] _? (All) 

f(a 0 ) - -f(ao) (rru 0 [y«ra 0 ♦ mKf/<T u ]- 1 (A12) 

where Kf is an emnirical fracture toughness oaraneter and m is a dimen- 
sionless emnirical coefficient which is not greater than unity. T'ote that 
if we let m«l and £7 U /7T/C n * ® Q * reduces to cq. (Al). 

G. Bo ckrath (1972 ) . - Equation (11) of ref. 5 for a center-crack rlat» 
ia limited to ^ D f 0,1 , which approximates an infinite plate. T hus 


, 2 

f ( a 0 ) * ^Tc~ («o)“ T**" 

(A13) 

f'(a G ) - -_Lf(, 0 )/a 0 

2*u> 

(All; ) 


(A1C 


where u> is a dimensionless empirical coefficient and K^ c is an empirical 
nnrameter having irrational units of (El, 00 ). 
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Crack Length, a 


a) Central representation of R-curve Instability ccnc*pt. 



Relative Crack Length, X 

b) R-curve Instability for a wide range of initial crack lengths. 


Figure 1. - R-curve Instability concepts. 
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Fracture Stress, £T 
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gure 5. - Actual R-curves and equivalent R-curvos (ERC) for hyocthetical material Urobtainiun 
infinite-width series. 
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Residual strength, equivalent R-curves, and estimated R-curve Tor 221 C -T?7 aluminum alloy sheet 



